On the sign of kurtosis near the QCD critical point 



M. A. Stephanov 
Department of Physics, University of Illinois, Chicago, Illinois 60607, USA 

We point out that the quartic cumulant (and kurtosis) of the order parameter fluctuations is uni- 
versally negative when the critical point is approached on the crossover side of the phase separation 
line. As a consequence, the kurtosis of a fluctuating observable, such as, e.g., proton multiplicity, 
may become smaller than the value given by independent Poisson statistics. We discuss implications 
for the Beam Energy Scan program at RHIC. 
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INTRODUCTION 

Mapping the QCD phase diagram as a function of 
temperature T and baryochemical potential /is is one 
of the fundamental goals of heavy-ion collision experi- 
ments. QCD critical point is a distinct singular feature 
of the phase diagram. It is a ubiquitous property of QCD 
models based on the chiral symmetry breaking dynamics 
(see, e.g., Ref.[l| for a review and further references). Lo- 
cating the point using first-principle lattice calculations 
is a formidable challenge (see, e.g., Ref.pj for a recent 
review and references). If the critical point is situated in 
the region accessible to heavy-ion collision experiments 
it can be discovered experimentally. The search for the 
critical point is planned at the Relativistic Heavy Ion 
Collider (RHIC) at BNL, the Super Proton Synchrotron 
(SPS) at CERN, the future Facility for Antiproton and 
Ion Research (FAIR) at GSI, and Nuclotron-based Ion 
Collider Facility (NICA) in Dubna (see, e.g., Ref.Q). 

The characteristic feature of a critical point is the di- 
vergence of the correlation length £ and of the magnitude 
of the fluctuations. The simplest measures of fluctuations 
in heavy-ion collisions are the variances of the event-by- 
event observables such as multiplicities or mean trans- 
verse momenta of particles. The singular, critical con- 
tribution to these variances diverges as (approximately) 
£ 2 , and would manifest in a non-monotonic dependence 
of such measures as the critical point is passed by during 
the beam energy scan 0, Q • L~l realistic heavy ion colli- 
sion the divergence of £ is cut-off by the effects of critical 
slowing down [{§, 6], and the estimates of the maximum 
correlation length are in the range of at most 2 — 3 fm, 
compared to the natural 0.5 — 1 fm away from the critical 
point. However, higher, non-Gaussian, moments of the 
fluctuations depend much more sensitively on £, accord- 
ing to Ref.p}. For example, the 4-th moment grows as 
£ 7 near the critical point, making it an attractive exper- 
imental tool. In this paper we follow up on the results 
of Ref. [7[ to point out that the sign of the 4-th moment 
could be negative as the critical point is approached from 
the crossover side of the QCD phase transition. 

The sign of various moments have been discussed in 
the literature in related contexts: see, e.g., discussion of 
the sign of the 3-rd moment in Ref. [H or the 6-th and 8-th 
moments in Ref. Q and also numerical lattice calculations 



in Ref. [1(3] where the possible sign change of kurtosis is 
noted. 

In this paper we shall address specifically the sign of 
the 4-th moment (or kurtosis) and do it in a more uni- 
versal and quantitative way than has been done previ- 
ously, by using the known parametric form of the univer- 
sal equation of state near the critical point. We empha- 
size universality of the behavior of the kurtosis and draw 
experimental consequences from these results. 



KURTOSIS AND UNIVERSAL EFFECTIVE 
POTENTIAL 

Let us begin, as in Ref. [l| , by describing fluctuations of 
the order parameter field er(x) near a critical point using 
the probability distribution 



P[*]~exp{-Sl[j]/T}.. 



(1) 



where f2 is the effective action (free energy) functional for 
the field a, which can be expanded in powers of a as well 
as in the gradients (we chose a — at the minimum): 
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Calculating 2-point correlator (a(x)a(0)) we find that 
the correlation length £ = m~ 1 . For the moments of the 
zero momentum mode ay = Jd 3 xa(x) in a system of 
volume V we find at tree level 
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where (a v ) c = 



— 3((iy) 2 denotes the connected 4-th 



central moment (the 4-th cumulant). The critical point 
is characterized by £ — > oo. The central observation in 
Ref.Q was that the higher moments (cumulants) K3 and 
K4 diverge with £ much faster than the quadratic moment 
«2- Here we shall point out that the sign of the 4-th 
moment K4 is negative in a certain sector near the critical 
point. More precisely, the 4-th cumulant is negative when 
the critical point is approached from the crossover side. 
Let us demonstrate this in several complementary ways. 

A simple way to see why the kurtosis is negative is by 
following the evolution of the probability distribution of 
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ay as we approach the critical point along the crossover 
line. In Ising scaling coordinates: along H = 0, t > ray. 
Away from the critical point, more precisely for £ 3 -C V , 
the central limit theorem dictates that the probability 
distribution of ay is Gaussian, with a vanishingly small 
kurtosis. As we approach the critical point the distri- 
bution develops non-Gaussian shape. This intermediate 
shape is a deformation of the Gaussian towards a two- 
peak distribution, corresponding to the phase coexistence 
on the opposite, first-order transition side (t < 0) of the 
critical point. Such a shape is clearly less "peaked" than 
the Gaussian, and thus corresponds to negative kurtosis. 

More quantitatively, the kurtosis vanishes as \/V at 
(almost) any point away from the critical point, i.e., 



K=K A /Ki=o{e/v). 



(4) 



The exception is the coexistence line (H — 0, t < ray). 
The distribution there has two peaks of equal height and 
its kurtosis is K = — 2+0(£ 3 /V). 1 It is important to note 
that this is only true strictly on the coexistence line H = 
0, for the moments measured around the symmetric point 
of the probability distribution of ay, which is actually 
a dip, not a peak, for t < 0. At any point close to the 
coexistence line, i.e., at H ^ 0, t < 0, the kurtosis around 
the dominant peak is positive. 

In the scaling regime (close to, but not at the critical 
point) where £ is much greater than the microscopic scale, 
a, but still much less than the linear size of the system: 
a < £ < V 1/3 , the coefficient of £ 3 /V in Eq. (gj can be 
expressed in terms of the couplings Aj using Eqs. ([3]): 
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(5) 



These couplings, and in fact the shape of the effective 
potential, is also universal. In particular, A4 scales with £ 
as A4 = where the universal value of A4 is known 

approximately to be 4.0 on the crossover line (see, e.g., 
Ref. llj for a review). 2 Since on the crossover line A3 = 
and A4 > 0, it is clear from Eq. © that K < 0. 

Away from the crossover line (H = 0, t > ray) the 
distribution is skewed: A3 ^ 0. This makes the kurtosis 
positive, according to Eq. (0, except for a certain sector 
around the crossover line. 



1 In fact, this transition of the shape of the distribution around a 
critical point is universal to all critical points of the same (Ising 
model) universality class. At finite V, at t = 0, the value of K 
is independent of V (the correlation length § ~ V 1 / 3 is as large 
as it can be at given volume). This value of K is a universal 
number (it depends only on the boundary conditions for a given 
universality class) and is well-known. It is usually expressed 
as the value of the Binder cumulant B4 1.6, which means 
K = B4. - 3 ss -1.4 < 0. 

2 As in Ref. , for simplicity and consistency with our overall level 
of precision, we neglect the anomalous scaling dimension rj, which 
is only of order few percent. 



THE UNIVERSAL EQUATION OF STATE 
AROUND THE CRITICAL POINT 

To extend this analysis away from the crossover line, 
i.e., to take into account A 3 ^ in Eq. ([5]), we need to 
know the equation of state, in particular, function 
of both Ising variables: reduced temperature t and mag- 
netic field H . In the scaling regime near t = H = this 
equation of state is also universal. For the Ising model it 
is known to order e 3 in the epsilon expansion as well as 
numerically. 

Before we discuss this universal form, let us keep in 
mind that the mapping of QCD phase diagram in the 
T, fig plane into t, H plane is not universal. However, 
this mapping is analytic, i.e., both functions t(T, hb) 
and H (T, /is) are analytic at the critical point, which is 
mapped into the origin, t(T c P,^ p ) = H(T cp ,lf^) = 0. 

The standard parametrization, Ref. [13], of the equa- 
tion of state in the scaling domain near the critical point 
is in terms of two new scaling variables R and 9 (it has 
been applied in the context of QCD before, Ref.[IH). De- 
noting the "magnetization" by M = (ay)/V we define R 
and 8 as 



m = R p e, t = R(i- e' z ), 



(6) 



Then the equation of state can be expressed in terms of 
the single function h(8) as 



H = R f3S h{9) 



(7) 



Unlike the explicit form of the singular equation of state 
M = M(t,H), the function h(6) is analytic. It has two 
zeros. One, at 6 = 0, corresponds to the crossover line 
(t > 0, H = ray), another, at some 8 = 6\ > 1, cor- 
responds to the coexistence (first-order transition) line 
(t < 0, H = ray). The function h{9) must also be odd 
since M(—H) = —M(H). The simplest function obeying 
all these requirements is a cubic polynomial 



h{6) =9(2 -3d 2 ). 



(8) 



where the value 6\ = a/3/2 is a good approximation to 
the universal value for the Ising model (and correct up to 
0(e 2 )). The choice ([8]) is known as the linear parametric 
model, Ref. 12]. It describes the equation of state with 
precision quite sufficient for our purposes. The linear 
parametric model is also known to be exact up to C(e 3 ). 

Using this parametric equation of state, we can cal- 
culate the moments by taking derivatives at fixed t, up 
to an overall normalization, unimportant in the present 
context (it can be fixed by Eq. (j3j). In particular, 



«4(t, H) 
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For our purposes, it would be sufficient to use the ap- 
proximate rational values of critical exponents (3 = 1/3 
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and 6 = 5, which are within few percent of their exact 
values in three dimensions. The result of Eq. ([9]) can then 
be simplified to 

, rrN 81 - 7836> 2 + 1056» 4 - 56» 6 + 2(9 8 , . 

n^H) = -12 ^3(3-^)3(3 + ^)5 • (10) 

We represent Ki{t,H) graphically as a density plot in 
Fig. [TJ We see that the 4-th cumulant (and kurtosis) 
is negative in the sector bounded by two curved rays 
H/t 133 = ztconst (corresponding to 6 w ±0.32). 




(a) 




(b) 



FIG. 1: (color online) (a) - the density plot of the function 
K4(t, H) given by Eq. (|10[) obtained using Eq. |9| for the linear 
parametric model Eqs. @, 0, © and /3 = 1/3, 5 = 5. The 
K4 < region is red, the K4 > - is blue, (b) - the dependence 
of K4 on t along the vertical dashed green line on the density 
plot above. This line is the simplest example of a possible 
mapping of the freezeout curve (see Fig. [2}. The units of t, 
H and K4 are arbitrary. 

Also in Fig. [1] we show the dependence of K4 along a 
line which could be thought of as representing a possible 
mapping of the freezeout trajectory (Fig. [2]) onto the tH 
plane. Although the absolute value of the peak in K4 
depends on the proximity of the freezeout curve to the 
critical point, the ratio of the maximum to minimum 
along such an H = const curve is a universal number, 
approximately equal to —28 from Eq. ([TU]) . 
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FIG. 2: A sketch of the phase diagram of QCD with the freeze- 
out curve and a possible mapping of the Ising coordinates t 
and H. 



The negative minimum is small relative to the positive 
peak, but given the large size of the latter, Ref.@, EH, 
the negative contribution to kurtosis may be significant. 
In addition, the mapping of the freezeout curve certainly 
need not be H = const, and the relative size of the posi- 
tive and negative peaks depends sensitively on that. 

The trend described above appears to show in the re- 
cent lattice data, Ref.[l(|, obtained using Pade resum- 
mation of the truncated Taylor expansion in (13- As the 
chemical potential is increased along the freezeout curve, 
the 4-th moment of the baryon number fluctuations be- 
gins to decrease, possibly turning negative, as the critical 
point is approached (see Fig. 2 in Ref.[10|). 

Another observation, which we shall return to at the 
end of the next section, is that — K4 grows as we approach 
the crossover line, corresponding to H = 0, t > on the 
diagram in Fig. 



1(a) 



On the QCD phase diagram the 
freezeout point will move in this direction if one reduces 
the size of the colliding nuclei or selects more peripheral 
collisions (the freezeout occurs earlier, i.e., at higher T, 
in a smaller system). 



EXPERIMENTAL OBSERVABLES 

In this section we wish to connect the results for the 
fluctuations of the order parameter field a to the fluctua- 
tions of the observable quantities. As an example we con- 
sider the fluctuations of the multiplicity of given charged 
particles, such as pions or protons. 

For completeness we shall briefly rederive the results of 
Ref.Q using a simple model of fluctuations. The model 
captures the most singular term in the contribution of the 
critical point to the fluctuation observables. Consider a 
given species of particle interacting with fluctuating crit- 
ical mode field ex. The infinitesimal change of the field Sa 
leads to a change of the effective mass of the particle by 
the amount 5m = gStr. This could be considered a def- 
inition of the coupling g. For example, the coupling of 
protons in the sigma model is gapp. The fluctuations Sf p 
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of the momentum space distribution function f p consist 
of the pure statistical fluctuations 6f p around the equi- 
librium distribution n p for a particle of a given mass, 
which itself fluctuates. This gives 

9n E 



6f * = 6f r + d£ 96a 



(11) 



Using this equation we can calculate the most singu- 
lar contribution from the critical fluctuations to the mo- 
ments or correlators of 5f p . The fluctuation of the mul- 
tiplicity N — V d f f p is given by 



SN = 5N Q + VgSad [ 



(12) 



where d is the degeneracy factor (e.g., number of spin or 
charge states of the particle) . Neglecting, for clarity and 
simplicity, the effects of quantum statistics, i.e., assum- 
ing n p <§C 1, we can use Poisson statistics for SN°. Using 
additivity of the cumulants (their defining property), and 
assuming SN a and 5a are uncorrelated, the contribution 
of the critical fluctuations can be expressed in terms of 
the corresponding moments of the critical field a fluctua- 
tions. For example, the contribution to the 4-th moment 
can be expressed as (cf. Refs.0, EH) 
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where 7 P = (dE p / dm)^ 1 is the relativistic gamma- factor 
of a particle with momentum p and mass m. The first 
term on the r.h.s. of Eq. (Ti~3"|) is the Poisson contribution. 
We neglected n p <C 1 in the quantum statistics factor 
(l±n p ) for simplicity, and we denoted by ". . ." other con- 
tributions, less singular at the critical point. The model 
is admittedly crude, but it illustrates the mechanism and 
correctly captures the most singular contribution near 
the critical point. 

In the region near the critical point where ka = {<Jy) c 
is negative, the 4-th cumulant of the fluctuations will 
be smaller than its Poisson value, (N). The measure 
defined in Ref.Q as lo 4 {N) = ((5iV) 4 ) c /(iV) will be less 
than 1. By how much will depend sensitively on the 
correlation length (as £ 7 ), i.e., on how close the freezeout 
occurs to the critical point, as well as on other factors 
(for protons, most significantly, on the value of /ib-) We 
shall not attempt to estimate this effect quantitatively in 
this paper. The analysis of Ref.[IH suggests, however, 
that this effect for protons can be significant compared 
to the Poisson value already for £ ~ 2 fm. 

Usual caveats apply: other (non-trivial) contributions 
to moments which do not behave singularly at the critical 
point can turn out to be relatively large. It is beyond the 
scope of the paper to estimate these effects. The size of 
these background contributions could, in principle, be 
determined experimentally by performing measurements 
away from the critical point. 



We conclude by asking an obvious question: has the 
effect of the negative kurtosis been observed? Data from 
STAR indicate that at y/s = 19.6 GeV the ratio k 4 /k 2 
might be substantially smaller than its Poisson value 1, 
see Fig. 6 in Ref . [l(| . while it is very close to 1 at higher 
yfs (smaller hb)- Unfortunately, the statistics accumu- 
lated in the short run at y/s — 19.6 GeV is clearly not 
sufficient to make a reliable conclusion. It would be in- 
teresting to see if this effect persists with more statistics 
at this energy. If confirmed, this result could indicate 
that the critical point is close, at somewhat larger values 
of [ib (smaller yjs). In this case, as we already discussed 
at the end of the previous section, the universality would 
also predict that the negative kurtosis effect should in- 
crease in more peripheral collisions at the same yfs. At 
smaller values of ^/s the effect should change sign, in- 
creasing kurtosis above its Poisson value. 
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